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ABSTRACT 


Thermal  flux  radiated  from  an  isothermal  slab  of  air  is  calcu¬ 
lated,  assuming  one -dimensional  radiation  flow  for  a  wide  range  of 
gas  temperatures  and  densities  and  using  the  latest  absorption  co¬ 
efficient  data  available. 

The  effects  of  assuming  an  optically  thin  g.^>n  and  different 
absorption  coefficient  data  are  considered.  Criteria  for  specify¬ 
ing  air  as  optically  thin  or  optically  thick  are  included.  The 
energy  lost  from  the  system  by  radiation  is  compar^:d  with  the  total 
energy  of  the  system. 

1.  INTRODUCTION 


Radiant  heat  transfer  can  be  an  Important  and,  in  some  cases, 
the  predominant  means  of  heating  a  missile  or  satellite  in  flight 
(ref  1).  Realization  of  this  effect  has  led  to  fundamental  studies 
of  the  radiative  properties  of  the  various  species  comprising  air 
(ref  2,  3).  Use  of  the  knowledge  obtained  from  these  studies,  how¬ 
ever,  is  restricted  by  the  complicated  nature  of  the  physical  sit¬ 
uation  and  the  resulting  mathematical  expressions.  Many  sinplify- 
ing  assumptions  are  required  before  a  calculation  can  be  made. 

This  report  Includes  calculations  of  the  energy  radiated  by 
air  over  a  wide  range  of  temperatures  and  densities.  The  deriva¬ 
tion  of  the  equation  used  in  the  calculation  is  outlined  explicitly, 
giving  the  assumptions  involved.  The  effects  of  using  different 
assumptions  are  shown.  Thus,  the  final  result  is  not  only  a  set 
of  data  based  on  the  most  recent  information  but  also  allows  an 
estimate  of  the  data's  validity  within  the  assuiuptlons  made. 

2.  RADIANT  ENERGY  TRANSFER  PROPERTIES 

2.1  Calculation  of  Thermal  Radiant  Energy 

The  radiation  emitted  by  high- temperature  air  is  calcula¬ 
ted,  starting  from  fundamental  principles  tliat  are  used  with  the 
most  recent  data  available  on  the  emissive  properties  of  air.  The 
independent  variable  is  the  temperature  with  density  ratio  as  a  par¬ 
ameter.  This  arrangement  allows  direct  entry  into  tables  of  shock- 
wave  properties,  such  as  those  calculated  in  references  3  and  4. 

The  radiant  flux  emitted  by  a  volume  of  high- temperature 
gas  is  determined  by  the  balance  of  emission  and  absorption  processes 
occurring  within  the  gas.  In  this  report  scattering  processes  are 
assumed  to  be  negligible.  Using  the  appropriate  expressions  defined 
in  appendix  A,  the  radiative  energy  change  in  a  volume  of  radiating 
gas  is  found  to  be 


5 


(1) 


dl 
_ u 

ds 


-  I  K  )  p 
V  u/  ^ 


(The  symbols  used  in  this  report  are  explained  in  appx  A,) 

and  are  dependent  upon  the  direction  of  the  radiation  flow. 

Assuming  local  thermodynamic  equilibrium,  we  can  write  for  the 
steady-state  case 

dl 

— p  =  B  -  I 
pKy  ds  u  i) 

where  the  effects  of  induced  emission  have  been  taken  into  account 
by  introducing 


K 


/ 
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K 
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(2) 


(3) 


Goulard  (ref  5)  and  Pai  (ref  6)  have  solved  equation  (2)  for 
a  one -dimensional  case.  Using  spherical  coordinates,  the  intensity 
is  then  independent  of  the  azimuth  angle  cp  but  dependent  on  the 
altitude  angle  G.  Using  their  results  (ref  5,  6),  together  with 
the  assumptions  of  an  isothermal  slab  of  a  gray  gas  (K  independent 
of  w)  with  no  interacting  sotid  surfaces,  we  find  for  l^he  emitted 
flux : 


Q  =  tT*  [1  -  2  E3(t)l  (4) 

where  t  is  the  integrated  optical  path. 

The  Planck  mean  absorption  coefficient  Kp  is  used  in  this  re¬ 
port.  It  is  the  mass  absorption  coefficient  K  weighted  with  re¬ 
spect  to  the  Planck  emissive  function.  (Appendix  A  defines  K  in 
detail.)  The  Planck  coefficient  is  independent  of  u  and  fulfflls 
the  requirement  of  a  gray  gas.  Kp  is  exact  only  for  an  optically 
thin  gas.  The  mass  absorption  coefficients  for  air  in  equilibrium 
have  been  calculated  by  Meyerott,  et  al  (i-ef  7),  using  quantum 
theoretical  and  experimental  techniques.  The  work  of  Gilmore  (ref  8) 
was  used  by  Meyerott  as  the  basis  for  determining  the  constituents 
of  air.  Gilmore's  constituent  results  do  not  differ  appreciably 
from  those  of  Hilsenrath  (ref  9,  10),  which  were  used  in  the  shock- 
wave  calculations  of  references  3  and  4. 

Armstrong  (ref  11),  using  the  work  of  Meyerott  (ref  7),  calcu¬ 
lated  Kp  as  a  function  of  density  and  temperature.  For  use  in  this 
report,  Armstrong's  tabular  results  have  been  carefully  plotted  in 
figure  1  to  give  a  continuous  set  of  values. 

The  function  B  (t)  appearing  in  equation  (4)  was  calculated 
by  Case  (ref  12).  For  use  in  this  report.  Case's  data  were  fitted 
to  the  curve 
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Planck  Mean  abaorption  coefficients. 


Eg(t)  =  a  +  bt  r  ct®  +  dt®  +  et* 


(5) 


Over  a  range  of  0.0 


<  t  ^  3.4, 


the  coefficients  were  found  to  be 


a  =  5.000000 
b  =-0.765638 
c  =-0.501736 
d  =-0.151243 
e  =  0.016886 


The  result  of  this  fit  is  shewn  in  figure  2  where  [1-2E  (t)j  is 
plotted  against  t  for  the  original  and  fitted  data. 

The  radiation  from  air  over  a  range  of  3000Oto  12,OOOPK  and 
densities  from  10~®  to  10  times  atmospheric  wore  calculated,  using 
the  mean  absorption  coefficients  of  Armstrong  (ref  11)  and  equa¬ 
tion  (4).  All  calculations  were  programed  by  the  author  for  the 
IBM  7090  computer.  The  results  are  shown  in  figure  3  for  a  path 
length  of  1  cm.  It  can  be  seen  that  the  radiation  emitted  from 
hot  air  can  be  appreciable,  even  at  the  lower  densities. 


2.2  Comparison  of  Calculated  Data 


The  enthalpy  H,  which  must  flow  into  a  gaseous  system  to 
reach  its  high-energy  state,  can  be  calculated  from 

H  =  pUgH'  (6) 

This  quantity  was  calculated,  using  the  de»*sity  and  shock  velocity 
found  in  reference  4,  and  the  enthalpy  per  unit  mass  of  reference  10. 

The  results  are  given  in  figure  3,  where  1  percent  of  the  inflowing 
enthalpy  is  shown.  It  can  Vi  seen  that  over  the  range  considered, 
the  radiated  energy  is  a  small  fraction  of  the  total  energy  of  the 
system.  Corrections  to  the  thermodynamic  properties  for  the  radiated 
enet^y  would  be  small. 

Figure  3  also  gives  the  energy  emitted  by  a  black  body  over 
the  temperature  range  considered.  It  can  be  seen  that  the  emitted 
energy  approaches  this  limit  at  the  higher  densities  and  temperatures. 

Kp  is  exact  for  an  optically  thin  case,  but  its  use  is  questionable  in 
the  high-density  and  high -temperature  regions.  An  estimate  of  how 
much  error  results  from  the  use  of  the  radiation  curve  based  on  K-  can 
be  made  by  examining  the  work  of  Strack  (ref  13),  who  evaluates  the 
integral 

Q  =  7?  B,.  L1-2E„(oK  s)j  dv  (7) 

J  t/  *>  V 

o 
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Figure  2.  Coapar 


using  th«  dsts  of  reference  7.  Strsck  finds  that  his  results  differ 
by  60  percent  from  those  obtained  uaing  Kp  in  the  region 


10-a  <  K  s  <  10  (8) 

-  P  - 

Strack's  results  are  compared  with  the  calculations  made  in 
this  report  (fig.  4)  for  densities  of  10  1,  and  10'”  times  atmospheric. 

There  is  no  difference  in  the  radiation  emitted  at  the  lowest  density, 
but  differences  up  to  35  percent  are  noted  at  the  higher  densities  in 
the  cases  considered. 

Strack's  results  are  effective  criteria  for  the  ranges  of 
optically  thin  or  optically  thick  air.  This  can  be  seen  by  consider¬ 
ing  figure  2.  For  t  >3.5,  the  quantity  [l-2l2(t)J  0.99.  Since  the 

function  approaches  1.0  asymptotically,  those  gases  emittiii99  percent 
of  the  black-body  value  itay  be  considered  as  black  bodies. 

The  optically  thin  gas  limit  can  be  obtained  from  equation 
(4),  using  the  expansion  found  in  reference  14: 

EgCt)  =  I  -  t  -  f(t)”  (9) 

which,  for  t  «  I,  yields 

Q  =  2t<ft^  (10) 

The  limits  obtained  above,  namely: 

10"”  >  K  s  >3.5  (11) 

-  p  - 

disagree  with  those  of  St rack  s  limits,  but  they  do  indicate  that  the 
differences  he  finds  are  in  the  region  where  the  gas  is  neither  opti¬ 
cally  thin  nor  optically  thick. 

The  criteria  for  optically  thin  or  optically  thick  gases, 
expressed  by  the  Inequality  given  in  equation  (11),  are  applied  to 
the  calculation  made  in  this  report  by  the  horisontal  lines  shown  in 
figure  1. 


The  effect  of  absorption  can  be  seen  by  comparing  the  cal¬ 
culations  already  made  with  the  optically  thin  case  as  given  by  equa¬ 
tion  (10).  Armstrong's  K  's  (ref  11)  were  used  in  equation  (10) 
and  the  results  shown  in  ^figure  4.  The  amount  of  radiation  is  the 
same  for  the  two  cases  at  the  lowest  density.  At  the  other  densities, 
it  can  be  seen  that  the  decrease  in  emitted  flux  caused  by  absorption 
is  about  10  percent  at  the  lower  temperatures.  The  flux  emitted  by 
the  optically  thin  gas  will  not  "cut  off"  at  the  black-bodv  limit  as 
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temperature  is  increased.  This  point  is  specifically  made  in  figure  4 
by  the  curve  representing  radiation  from  optically  thin  gas  at  the 
highest  density. 

The  radiant  flux  from  an  optically  thin  gas  has  been  calcu¬ 
lated  by  Kivel  (ref  1),  who  uses  his  own  calculation  (ref  15)  for 
the  emissive  properties  of  air.  His  results  are  compared  with  the 
optically  thin  calculation  made  in  this  report  in  figure  4.  The  lack 
of  agreement  shown  is  chiefly  due  to  the  choice  of  different  emissive 
properties  of  air  for  each  calculation. 

3.  SUMMARY  AND  CONCLUSIONS 


All  the  calculations  siade  in  this  report  assume  an  isothermal 
slab  of  gas  that  is  in  a  state  of  local  thensodynamic  equilibrium.  It 
is  further  assumed  that  no  solid  surfaces  are  interacting  with  the 
radiating  gas  and  that  tlie  radiation  flow  is  one-dimensional. 

Within  the  limits  imposed  by  these  assumptions,  equation  (4)  is 
the  most  general  expression  used.  This  equation  is  used  to  calculate 
the  radiation  from  a  slab  of  air  1.0  cm  thick  over  a  wide  ranf^e  of 
temperatures  and  densities  (fig.  3).  Radiation  from  slabs  of  differ¬ 
ent  thickness  may  be  calculated  using  the  data  of  figure  2. 

The  results  of  using  equation  (4)  with  K  are  not  completely  ac¬ 
curate  because  K  is  exact  only  in  optically  ?hin  gases.  Corrections 
may  be  introduceS  by  using  data  derived  from  equation  (7).  These  cor¬ 
rections  range  up  to  60  percent  of  the  total  radiation  but  apply  only 
over  limited  ranges  of  temperature  and  density.  A  conservative  esti¬ 
mate  of  the  regions  where  these  corrections  apply  is  given  by  the  in¬ 
equality  expressed  in  equation  (8).  This  inequality,  together  with 
the  data  of  figure  1,  can  also  be  used  to  define  regions  where  air  is 
optically  thin  or  optically  thick. 

The  enthalpy  represented  by  the  escaping  radiation  is  compared 
with  the  total  enthalpy  of  the  system  (fig.  3).  Within  the  tempera¬ 
ture  range  considered,  the  escaping  radiation  has  little  effect  or 
the  energy  state  of  tiie  system. 

The  assumption  of  an  optically  thin  gas  is  often  made.  Equation 
(10)  is  the  result  of  this  assumption.  This  result  is  compared  with 
the  case  of  appreciable  self-absorption  as  represented  by  equation  (4) 
in  figure  4.  A  difference  of  approximately  10-percent  exists  between 
the  two  cases  at  lower  temperatures.  The  optically  thin  expression 
does  not  "cut  off"  at  the  black-body  limit,  so  the  divergence  between 
the  two  cases  steadily  Increases  as  temperature  increases. 

Although  the  absorption  coefficient  used  is  based  on  the  most 
recent  data  available  (ref  11),  a  comparison  with  older  but  accepted 
data  (ref  1)  indicates  that  more  research  is  needed  on  the  fundamental 
properties  of  air. 
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APPENDIX  A.  NOMENCLATURE  AND  DEFINITIONS  USED  IN  CALCULATION 


OF  RADIANT  ENERGY  TRANSFER 


Planck’s  Function,  By 


=  (2h/c2)  \p 


exp 


(A-1) 


equals  the  energy  radiated  from  a  point  in  a  black  body  in  a  given 
direction  per  unit  solid  angle  per  unit  time  per  unit  projected 
area  of  the  emitter. 


c  =  velocity  of  light, 
h  =  Planck’s  constant, 
k  =  Boltzmann’s  constant, 

T  =  absolute  temperature,  and 

U  =  frequency  of  emitted  radiation;  using  u  as  a  subscript 
denotes  a  frequency  dependent  quantity. 

B  =  oT^  (A-2) 

equals  energy  radi"  ted  from  a  point  in  a  black  body  at  all  frequen¬ 
cies  into  a  solid  angle  of  2v  per  unit  time  per  unit  projected  area 
of  the  emitter.  a  is  the  Stefan-Boltzmamconstant . 

OD 

E  (t)  =  j  u"®  [exp  (-tu)]du 
o  J 

i 

u  =  sec  0  (A-3) 


Specific  Intensity,  ly 


I 

u 


lira 

(dj,  du),  dt,  du 


0) 


dE 

_ u _ 

dcT  cos  0  duj  dt  du 


(A-4) 


dE  =  total  amount  of  energy  passing  through  tne  element  of 
^  area  da  inside  the  element  of  solid  angle  duj  in  time 

interval  dt  and  in  the  frequency  interval  (u,  u  +  du). 

0  =  angle  between  the  normal  to  the  emitter’s  surface  and 

the  given  direction. 

Mass  (or  Monochromatic)  Absorption  Coefficient,  Ky 


dl 


u,ab 


-K  0  I  ds 
u  u 


(A-5) 
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dl  ^  =  '-adiative  intensity  absorbed  in  gas 
w^ab 

ds  =  element  of  path  length  over  which  radiation  travels. 
0  =  gas  density. 


Mass  Absorption  Coefficient  Considering  Induced  Emission,  Ky 


K'  =  K  1-exp 
u  u 


- 

kT/. 


(A-6) 


Emission  Coefficient,  jy 


dE  .  =  i  dm  du)  dt 

i;,-tm 

energy  emitted  from  an  element  of  mass  dm  into  the 
element  solid  angle  du)  in  time  inten'al  dt  and  in 
the  frequency  interval  (u,  v  +  dy). 

K  B  (gas  in  thermodynamic  equilibrium) 

V'J 


(A-7) 


Flanck  Mean  Absorption  Coefficient,  K 


K  = 
P 


•no  k'  1  av 
^  Jo  u  V 

jT* 


(A-8) 


Flux  of  Radiation,  Q 


Q  =  r  i  I  dv  cos  0 

J  J  u 


(JU  o 

equals  net  flux  of  radiation  per  unit  area  per  unit  time 


(A-9) 


Q  H  €  B 
U  1’  u 


(A-10) 


€  =  emissivitv 

V 

Optical  path,  t 


■  r- 


pdsdy 


(A-11) 


Gas  Properties 


H  =  enthalpy  flux 

=  enthalpy  per  unit  mass  of  gas 
U  =  vslocilv  of 
p  =  gas  density 

Q  =  gas  density  at  1-atm  pressvire 


